S>-Ri37  129 

COUPLING  BETWEEN  TWO 
THROUGH  NULTILRVERED 

ARBITRARILV 

SHIELDS(U) 

ORIENTED  DIPOLES 
ILLINOIS  UNI V  AT  URBANA 

ELECTROMAGNETICS  LAB 

R  VANG  ET 

AL  JUN  82  UIEM-82-7 

UNCLASSIFIED 

N00014-81-K-0245 

F/G  20/3 

NL 

MICROCOPY  RESOLUTION  TEST  CHART 
NATIONAL  BUREAU  OF  STANDARDS-1963-A 


ELECTROMAGNETICS  LABORATORY 
TECHNICAL  REPORT  NO.  83-7 


June  1983 


AD  Al  371  29 


COUPLING  BETWEEN  TWO  ARBITRARILY  ORIENTED 
DIPOLES  THROUGH  MULTILAYERED  SHIELDS 


R.  YANG 
R.  MITTRA 


TH~  r- ■' ,  .. 

TCC  "  •'  ■ 


uaval  , 

- - 


mic 

_ _ V,  JAM  2  31984 

This  document  has  been  approved  I  "  ;  - 

for  public  release  and  sale;  its  f 

distribution  is  unlimited. _ I  V 


ELECTROMAGNETICS  LABORATORY 
DEPARTMENT  OF  ELECTRICAL  ENGINEERING 
ENGINEERING  EXPERIMENT  STATION 
UNIVERSITY  OT^EBSIB**  URBANA-CHAMPAIGN 
URBANA,  ILLINGfCBttftJbfc, 


SUPPORTED  BY 
CONTRACT  NO.  N00014-81-K-0245 
PHYSICAL  SCIENCES  DIVISION 
OFFICE  OF  NAVAL  RESEARCH 
DEPARTMENT  OF  THE  NAVY 
ARLINGTON,  VIRGINIA  22217 


SKCUNtTV  CLASSIFICATION  OF  THIS  FAGC  'Wimi  Data  Sntarad; 


REPORT  DOCUMENTATION  PAGE 


RKFQRT  NUMSI 


OAfiC  READ  INSTRUCTIONS 

I  BEFORE  COMPLETING  FORM 

t  GOVT  ACCESSION  NOJ  3.  ACCIFlCNT'S  CATALOG  NUMSCR 


4.  TITLE  (m*  SuMMla) 

COUPLING  BETWEEN  TWO  ARBITRARILY  ORIENTED  DIPOLES 
THROUGH  MULTILAYERED  SHIELDS 


17.  AUTHOR^ 


R.  Yang  and  R.  Mittra 


5.  type  op  report  *  perioo  covereo 

Technical 

li  PERFORMING  ORG.  REPORT  NuMSER 

EM  83-7;  UILU-ENG-83-2546 

§.  contract  or  grant  NUMae^s) 

N00014-81-K-0245 


9.  PERFORMING  ORGANIZATION  NAME  AND  A00RES3 

Electromagnetics  Laboratory 

Department  of  Electrical  Engineering 

University  of  Illinois.  Urbana.  Illinois  61801 _ 

It.  CONTROLLING  OFFICE  NAME  AND  ADDRESS 

Office  of  Naval  Research 
Department  of  the  Navy 
Arlington,  Virginia  22217 

t*.  MONITORING  AGENCY  NAME  A  AOORSSSf//  dUtarantlrom  Controlling  Ottlca) 


to.  program  element,  project,  task 

AREA  A  WORK  UNIT  NUMBERS 

Project  No.  410 


12.  REPORT  DATE 

June  1983 

13.  NUMBER  OP  PAGES 

40 

15.  SECURITY  CLASS,  (ot  tht a  rmport) 


Unclassified 


ts«.  DECLASSIFICATION  DOWNGRADING 
SCHEDULE 


|  1ft.  DISTRIBUTION  STATEMENT  foi  thia  Raport) 


Distribution  Unlimited 


t7.  DISTRIBUTION  STATEMENT  'o  t  thm  aba  tract  ant  arid  in  Sloe*  20,  it  ditlarant  from  Raport) 


I  18.  supplementary  NOTES 


M9.  KEY  WORDS  on  ravaraa  a  Ida  it  nacaaaary  and  idantify  by  block  nunbtr) 


Multilayered  Shields;  Plane-Wave  Spectral  Representation;  Coupling; 
Arbitrarily  Oriented  Dipoles;  Transformed  Matrices 


20. ^ABSTRACT  fContinua  on  raaaraa  at  da  it  nacaaaary  and  idantify  by  block  numtiar) 

An  analysis  of  the  electromagnetic  coupling  between  two  arbitrarily  oriented 
dipoles  through  planar  multilayered  shields  is  made.  This  method  of  approach 
is  based  on  the  plane-wave  spectral  representation  of  radiation  fields  and 
the  use  of  transformation  matrices  for  multilayered  media.  Numerical  results 
are  found  to  be  in  good  agreement  with  the  experimental  data  obtained  in  the 
frequency  range  from  10  kHz  to  1  MHz.  y 


do  1473 


EOt T* ON  OP  ’  NOV  *5  %  OBSOLETE 


UNCLASSIFIED 


SSCjRl^v  CLASSlP  CATiQn  2 c  ”HiS  PAGE  .  Jmta  zntara*i) 


UILU-ENG-83-2546 


Electromagnetics  Laboratory  Report  No.  83-7 


COUPLING  BETWEEN  TWO  ARBITRARILY  ORIENTED  DIPOLES  THROUGH  MULTILAYERED  SHIELDS 


by 

R.  Yang* 
R.  Mlttra 


Electromagnetics  Laboratory 
Department  of  Electrical  Engineering 
University  of  Illinois  at  Urbana-Champaign 
Urbana,  Illinois  61801 

*0n  leave  from 
Xian  Jiaolong  University 
Xian,  Shaanxi  Province 
The  People's  Republic  of  China 


Technical  Report 


June  1983 


Supported  by 

Contract  No.  N000I4-81-K-0245 
Physical  Sciences  Division 
Office  of  Naval  Research 
Department  of  the  Navy 
Arlington,  Virginia  22217 


ABSTRACT 


An  analysis  of  the  electromagnetic  coupling  between  two  arbitrarily 
oriented  dipoles  through  planar  multilayered  shields  is  made.  This  method  of 
approach  is  based  on  the  plane-wave  spectral  representation  of  radiation  fields [1] 
and  the  use  of  transformation  matrices  for  multilayered  media.  Numerical 
results  are  found  to  be  in  good  agreement  with  the  experimental  data  obtained  in 
the  frequency  range  from  10  kHz  to  1  MHz. 
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1 .  INTRODUCTION 


It  is  well-known  that  an  arbitrary  field  can  be  represented  in  terms  of  its 
plane-wave  spectral  representation.  In  this  paper,  we  employ  this  represen¬ 
tation  to  express  the  incident  and  scattered  fields  from  multilayered  media  as 
follows 

E(x,y,z)  -  //f(a,6)  eJ<  0“+6,+Yz)  dads  (1) 

H(x,y,z)  -  //  »«,»  dadS  (2) 

where  y  m  A*  -J-  t  m  and  ~  on  top  indicates  Fourier  transform.  For 

a  plane  wave  that  is  incident  on  the  boundary  between  two  media  from  an 
arbitrary  direction,  it  is  possible  to  find  a  matrix  relationship  between  the 
fields  in  the  two  regions  separated  by  the  boundary.  The  matrix  may  be  termed 
the  transformation  matrix  and  will  be  very  useful  for  our  analysis.  In  order  to 
simplify  the  representation  of  the  transformation  matrix,  it  is  useful  to 
resolve  the  incident  wave  into  two  components,  one  parallel  to  the  plane  of 
incidence,  and  the  other  perpendicular  to  the  same.  These  two  components  will 
be  denoted  by  subscripts  (I  -  parallel,  1  -  perpendicular),  propagate  indepen¬ 
dently,  and  remain  uncoupled  throughout  the  processes  of  reflection  and 
transmission  at  the  interface  between  dissimilar  media.  Hence,  the  transfor¬ 
mation  matrix  for  the  fields  expressed  in  terms  of  the  parallel  and  perpen¬ 
dicular  components  (to  the  plane  of  incidence)  becomes  diagonal. 


2.  TRANSFORMATION  MATRIX  (Tj  FOR  PLANE  WAVES 


■M 

& 


In  this  section,  we  first  derive  the  transformation  matrix  for  a  single 
interface  between  two  dissimilar  media.  The  extension  to  the  multilayered  case 
is  considered  later.  Let  us  define  two  reference  planes  that  are  parallel  to 
the  boundary  at  Z  *  0,  namely,  plane  1  and  plane  2,  as  shown  in  Fig.  1.  The 

l®  l©  20  2© 

waves  E^  and  E^  are  at  plane  1,  and  E^  and  E^  are  at  plane  2.  The  ©  sign 
indicates  waves  propagating  in  the  +z-direction,  whereas  the  0  sign  indicates 
waves  propagating  in  the  opposite  direction. 
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Next  we  resolve  the  incident  electric  field  of  the  incident  wave  into  the 
parallel  and  perpendicular  components.  Figures  la  and  lb  illustrate  the 
geometry  of  propagation  of  the  perpendicular  and  parallel  components ,  respec¬ 


tively.  The  boundary  condition  that  must  be  satisfied  on  the  interface  between 
the  two  dissimilar  media  is  that  the  tangential  components  of  the  fields  E  and  H 
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be  continuous.  Enforcing  these  conditions  we  get 

Ef  eJrldl  -  Ef  -  Ef  -  Ef 

-l“  .Jr‘d‘  co.  6.  -  a®  e'Jrldl  CO,  9.  -  ,'Jr2d2  co.  92  *  af 


'1  1 


1  1 


cos  9„ 


where  dj  and  d 2  are  distances  of  the  terminal  planes  from  the  interface,  and 
0j  and  02  are  the  angles  between  the  directions  of  propagation  and  the  normals 
of  the  interface,  as  shown  in  Fig.  1.  Therefore,  we  have 


cos  0 


1  2  12 

tj  -  a  -  6 


1  Ik, 


COS  0„ 


-  a  -  6 


m 


*  -a.  ^  A  *  A.  ^  A  i  m m  _  •  _  " 


Note  that  the  phase  factor  e-1  is  assumed  here.  Rewriting  Eq.  (4)  and  using 
H  ■  /e/u  E,  we  obtain 


jr,d 


1“1  _ 


Ef.-3rldl) 


cos  9 
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j 

e 


r2d2 


)  COS  0 


(7) 


Multiplying  Eq*  (3)  by  (V^/u^  cos  )  and  subtracting  from  Eq, 


( 7 ) ,  we  have 


Similarly,  multiplying  Eq.  (3)  by  (/e^/ jj^  cos  9^ )  and  subtracting  from  Eq.  (7), 
we  have 


From  Eqs.  (8)  and  (9),  we  obtain 
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Ei 

E 

-  [T] 

Ei 

0le 

i 

29 

Ei 

V—  -J 

Ei 

(10) 


where 


(13) 


E  E 

These  matrices  [T]  and  [T]  are  used  to  calculate  the  coupling  between  two 
1  1 

magnetic  dipoles.  For  the  coupling  between  two  electric  dipoles,  it  is 

necessary  to  find  the  transformation  matrices  for  the  magnetic  field,  namely, 

H  H 

[TJ  and  t T]  •  It  is  evident  that  Fig.  1(a)  also  describes  the  propagation  for 
1  1 

the  parallel  components  of  the  magnetic  fields,  and  Fig,  L(b)  the  perpendicular 


components  of  the  magnetic  fields.  With  Eq.  (4)  and  E  replaced  by  H  in  Eq.  (3), 
It  is  not  difficult  to  show  that  the  transformation  matrices  for  magnetic 
fields  are  given  by 


where 


V  S 


Having  derived  the  matrices  for  a  single  interface  between  t  dissimilar 
media,  we  now  proceed  to  generalize  the  analysis  for  the  case  of  multilayered 


media*  Suppose  we  have  n  different  media  and  (n  -  1)  boundaries  as  shown  in 
Fig.  2.  From  the  preceding  results,  we  have 
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where  W  stands  for  the  waves  propagating  in  the  +z-direction  in  the 

t  h  jjQ 

i  medium,  and  W  stands  for  the  waves  propagating  in  the  -z-direction  in  the 
iC^  medium,  i.e. , 


W 


i® .  Ei®_  E»  H»  ot  Hi® 


wie  ■  if,  Ef ,  Hf .  or  Hf 
We  obtain  the  following  results: 
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n© 

n© 


i  ■  1,  2,  3,  ...  n 
i  ■  1,  2,  3,  ...  n 


8 
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(18) 


The  matrix  [ T J ,  can  be  termed  the  transformation  matrix  for  the  N-layer 
medium. 

It  is  obvious  that  the  transformation  matrices  only  depend  on  the  direc¬ 
tions  of  the  waves,  the  constitutive  parameters  of  the  media,  and  the  locations 
of  the  terminal  planes*  For  the  lossy  medium,  all  of  the  above  formulas  are 
still  applicable,  provided  we  replace  e  by  ( ef  +  je")* 


3.  SPECTRAL  REPRESENTATION  OF  THE  RADIATION  FIELDS 


The  fields  due  to  a  z-oriented  electric  dipole  located  at  (x'.y'.z')  can  be 
derived  from  a  z-oriented  electric  vector  potential  A^,  which  satisfies  the 
equation 


2  2 
7  A  +  It  A 
z  z 


-Idi  5(x  -  x')  <5(y  -  y')  6(z  -  z') 


(20) 


where  Idl  is  the  dipole  moment  of  the  electric  dipole.  Equation  (20)  can  be 
solved  via  the  Fourier  transform  technique  to  yield  (see  Appendix  I) 


A  (x.y.z) 
z 


Idi 

8„1 2 


j  (  a( x-x '  )+B(y-y'  )+y(z-z' )  ] 


dadg 


(21) 


-  B 


From  Maxwell's  equation  we  have 


r  3A 

Hx  "  17" 


H 


v  y 


1  g2 2 

1  z 


jue  3z  3x 


.  32A 

1 _ z 

jtue  3z3y 


(22) 


C‘ 


1  s2a* 
-j  OJUA  +  — - 1 

2  3z2 


Therefore,  the  various  field  components  for  a  z-oriented  electric  dipole  are 
given  as  follows 


H  »  rr  P 

X  8w2  11  A1  -a1-  62 


gj  [  a(x-x' )+6(y-y' )+y|z-z'  |]  dQdg 


„  .iil  ff  j[a(x-x')+B(y-y'HYl*~z*  |1  dad6 

y  8,2-i#T7T7 

E  -  Id*  r7z®  0(x~x,)+8(5^~y'  )+y|z-z  ’ll  dodg 

x  s  X2--^ 

E  ml*l  [IzA  gjtoCx-x^+eCy-y'J+Yl*-*'  |]  dadg 
y  8  it2  -i  “e 

2  2 

.  -iii  r”? .jlaC«-x'H6(y-y')+ir|s-t’ II 

-*  8.2  £~JZ  g2 


where 


-o-8 


Similarly,  we  can  derive  the  various  field  components  for  an  X-oriented 
electric  dipole  as  follows 


'H  -Ii|  //  ^Jlo(JMi,)+B(y-y*)+Yl*-«,|J  dad0 


8ir  -“ 


II 


8/  --  #  .  .  „2 

- j .  *  "(k2  “  a2) 

—  II - ===== 

8ff2  "•  aieA?  -  a2  -  62 


JIoCx-x'J+SCy-yM+Ylz-z’  |]  dadg 


JlaCx-x^+BCy-y'J+Ylz-z' I)  dadg 


gj [ o(x-x' )+8(y-y' )+y\z-z'  I )  dade 


8*  -  mA1  -  a2  -  e2 


.E  -  Id*  rrZ2  eJlaCx-x'HB^-y'J+Ylz-z' 

z  Sir2  --  “e 


dad8 
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r  / 


From  duality,  we  know  that 


(25) 


E  ♦  H  J  Id*  +  kdi  J  e  +  u 

H  ♦  -E  1  kdJt  ♦  Id*  1  u  ♦  e 

where  kd*  is  the  dipole  moment  of  the  magnetic  dipole.  Consequently,  the 
various  field  components  for  the  magnetic  dipole  can  be  obtained  directly  from 
those  of  the  electric  dipole  upon  application  of  the  duality  principle. 

The  various  field  components  for  the  z-oriented  magnetic  dipole  are 


r 


kd* 


If 


-$ 


Sir2  “•  A2  -  a2  ~  02 


ej  [  a(x-x*  )+0(y-y' )+y|z-z'  |]  dadg 


v  m  kd*  ff _ 2 _ (  o(x-x'  )+f3(y-y'  )+y  |z-z'  |  ] 

?  7?  H  /2  :  a2  _  2 


dad  B 


X  8ir  -•  *"** 


(26) 


H  -Mi  flzi  ejta(x-x,)+6(y-y,)+Y|z-z*  |] 

y  Sr2  —  “W 


2  ^  „2 


H  -  kd*  ff  ~(a  =  LJ  -l[a(x-x,H6(y-y,)+Yl2-2l  |1  dad 


Likewise,  the  components  for  the  x-oriented  magnetic  dipole  are  given  by 
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4.  NUMERICAL  COMPUTATION 


4.1.  Collinear  Magnetic  Dipoles  (or  Small  Coaxial  Loops) 

Suppose  we  have  two  z-oriented  collinear  magnetic  dipoles  separated  by  a 
shield.  Loop  1  is  a  transmitting  antenna  placed  at  the  origin  of  Cartesian 
coordinates  (0.0,0),  whereas  Loop  2  is  a  receiving  antenna  at  the  point 
(0,0,Zq),  as  shown  in  Figure  3. 

Letting  the  terminal  planes  at  z  *  0  and  z  ■  z,  respectively,  the  incident 
electric  fields  E^  ,  E^  for  the  magnetic  dipole  at  the  plane  1  are 

>(x,y,0)  -M;<  -t  dad6 

8*  A2  -  A  -  s2 


II - 2 - ( oort-ey)  iaii 

8"  A2-  A-  B2 


(28) 


From  the  matrix  Eq.  (A. 12)  in  Appendix  2,  we  obtain 


r , 

a 

~1© 

E1 

- 

A2  ♦  B2 

/a2  +  82 

E1® 

X 

pi© 

B1 

Jj  -  «2  -  B2 
,  /2  2 

bA2  -  a2  -  B2 
s  ,  /  2  2 

El® 

y 

W.  — 

^  k/a  +  8 

k/a  +  6 

(29) 


and 


E^x.y.O) 

l 


ltd  a 
8tt2 


Jl 

-V  a 


"/r 


±  6  j(ax+6y) 


dodS 


-  8 


(30) 


^E*®(x,y,0)  -  0 

This  result  shows  that  the  direction  of  every  incident  electric  field  at  plane  1 
is  always  perpendicular  to  its  plane  of  Incidence,  because  the  electric  field  of 


15 


the  z-oriented  magnetic  dipole  in  the  cylindrical  coordinates  (p,<J>,z)  has  the 
one  ^component  • 


w1®  -  e‘® 

W1®  .  E>® 


Finally, 


re 

W  -  E 


~1©  1  ~1© 

E  * -  E 

1  T1  1 

11 


EL  ER 


where  T  is  the  first  element  of  the  matrix  [T]  which  is  [T]  *[T]  .  The 

1 1  EL  1  11 

matrix  [T]  is  the  transformation  matrix  for  tne  left  boundary,  whereas 
ER  1 

the  [T]^  is  the  one  for  the  right. 

Therefore,  one  may  obtain  the  field  at  plane  2 


-2®,  v  fed  1  r  r 

Ei  (x,y.*0>  ■  — f  // 

8  it 


ej(«cey)  dod6 


-  a  -  0 


In  order  to  find  the  magnetic  fields  received  by  loop  2,  it  is  necessary  to 

OA  nm  vts 

transform  the  E  into  the  components  E  and  E  .  From  the  matrix  Eq.  (A. 13) 

1  x  y 

in  Appendix  2,  we  obtain 


(A13) 


E?(x>y’V  11 

x  u  8n  — 


j(ax+8y) 

e  dcd  8 


TjjA2  -  -  82 


<*.y.  V  ’  ^  // 

7  8  TT  -oo 


j(ooc+8y) 
e  dodg 


tiX?T7T7 


It  is  obvious  that  the  z-oriented  magnetic  dipole  placed  at  the  z-axis  can 
only  couple  the  z-component  of  the  magnetic  field.  From  Maxwell* s  equation  we 


H  -  j 
z  J 


(DU  3x  3y 

'w  w/ 


H  (x,y,i0)  -  — J—  // - — - — - eJ( coc^Sy)  date 

8^u  ^jsrrrr? 

In  order  to  evaluate  this  Integral,  let  us  make  the  following 
substitutions  [2]: 


o  ■  X  cos  ? 
8  ■  X  sin  5 


x  ■  p  cos 


y  ■  p  sin  <j) 


It  is  worth  noting  that  the  first  element  TT"  of  the  matrix  [T]  is  only  a  func- 

11  1 

tlon  of  X  and  independent  of  £,  because  the  factors  contained  in  the  matrix 


become 


cos  9 


2  -  a2  -  82  A2  -  I2 


c2  -  a2  -  82  -  A2  ~  X2 


Hence,  one  can  separate  the  above  double  integral  and  obtain 
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H  (  p»  z  ) 
z  0 


1  "  -X  dX  r  i  oX  cos  ( £-*) 

2 —  i  i  e  d5 


2  J 
8tt  ujy  0 


•f.^2  -  »2 


8"  * 0  Tf,x?r? 


JQ(pX)  dX 


Making  another  substitution 


X  »  sin  ? 


we  have 


.  .  kdi  ,  -2irk^  sin^  x,  r  .  .  , 

H_(  P»  ♦»*q)  "  5 / - j -  *  J0(kP  8in  0 

z  u  8vzwu  c  tf. 


where  c  is  a  contour  in  the  complex  5-plane,  as  specified  by  the  transformation 
X  *  k  sin  g.  In  this  case,  the  magnetic  dipole  is  placed  in  the  free  space,  so 
k  is  real.  Let  k  *  k'.  Then 


X  »  X'  +  jX"  *  k' (sin  ch?”  +  j  cos  she’ 


It  is  easy  to  show  that  the  mapping  of  Eq.  (34)  transforms  the  quadrants  of  the 
X-plane  into  parallel  strips  of  width  ir/2  radians.  The  path  of  the  integrands 
in  Eq.  (39),  that  runs  from  X  ■  0  to  X  ■  +«  in  the  right-half  plane,  is  now 
transformed  into  the  path  c,  as  shown  in  Fig.  4. 


Finally,  we  have 


V>>**o>  •yk<*ln  '■> d'’ 


Tj L ( cos  ?’) 


<■00  3 

_  f  slsi—  J0(kp  chC”)  d?-) 

00  x 


Tu(shC") 


The  magnetic  field  received  from  the  z-oriented  magnetic  dipole  placed 


along  the  z-axis  is 


Hz(0,0,z0) 


-k^kdJl 


/  Sin  Hr'  Ch  dg 


0  ml 


t:.(cos  ?•) 


T^Csh?") 


If  the  shield  is  removed,  it  becomes  a  free  space.  From  Eqs.  (26),  the  com¬ 


ponent  of  the  magnetic  field  on  plane  2  is  obtained 


H  .  JSji-  Z0  dads  (45) 


8ir  o)M  -oo  ^2  _  2  _  &2 


Using  the  same  substitutions  in  Eqs.  (37)  and  (40),  one  may  obtain 


-TTSST  >  slah'  Vk»  sl“  t’>  eJ"Z°C°8  d,:, 


+  /"  chV  J  (kpch?-)  ekZ°ShC"  d?’ 
0  U 


„  -k3kdJl  ,  3  ,  jkzncos  r' 

HZ  (0,0,.)  --47=T  /  s1"3  <  '  0  d<' 


/  chJ  5-  e 
0 


3 ...  kz0shc" 


v  V  «-  V  V  *v  ".\V/  \  V 

-  •  .  *  .  ■  *  «/  .  -  -  v  v  v  \  \  •  .  -  • 

*  •*,  *',  •*«**«/  */  ■/  *,*  /•  *  ■*  v  -«*  V- 


•*,  v.  a,  - 


Let  us  define  the  shielding  coefficient  of  the  shield  as 


lHzo(°»0>zn) | (free  space) 

S  -  20  log -  (db)  (48) 

|H  (0,0, z  )|  (with  shield) 
z  U 

It  can  be  used  to  describe  the  coupling  between  two  coaxial  z-oriented  small 
loops  separated  by  the  shield.  The  large  S  means  that  the  coupling  is  weak. 
4.2.  Collinear  Electric  Dipole 

Consider  two  z-oriented  collinear  electric  dipoles  separated  by  a  shield, 
as  shown  in  Fig.  5.  Invoking  duality  (Eq.  (25)),  it  is  possible  to  find  the  z- 
component  of  the  electric  field  at  plane  2  directly 


Ez(p,4>,z0) 


/  - ^ -  •  Jn(kpsin  5')  d z’ 

0  1  u 
T^cos  ?') 

-CO  3 

-  /  Ch  JQ(kp  chc")  dC") 

*  T^she”) 


-k2Idg 

2 

4  TT  UI€ 


(49) 


and 


E  (0,0, z  ) 
z  0 


k3Idf 

7~2 

4  ir  we 


I  sin3  dr1 

°  Tn(cos 


d{- 

T^(she") 


(50) 


But  t|.  is  the  first  element  of  the  transformation  matrix  for  the  perpendicular 
11  H  E 

component  of  the  magnetic  fields  [ T 1  in  Eq.  (15),  Instead  of  [T]  . 

1  1 


For  free  space,  we  have 


Similarly,  we  can  define  the  shielding  coefficient  of  the  shield  as 


S 


E 


IEzqCO.O.Zq) |(free  space) 
^  |Ez(0,0,z0) |(with  shield) 


(dB) 


(5. 


It  describes  the  coupling  between  two  collinear  electric  dipoles  separated  by 


shield 


5.  NUMERICAL  RESULTS 


We  have  two  coaxial  loops  separated  by  a  double  copper  shield.  Let  us  com¬ 
pute  the  shielding  coefficient  of  this  shield.  The  thickness  of  each  shield  is 
(1/1000)”;  the  separation  between  the  two  shields  is  ds  «  (1/2)"  or  (1/8)”;  and 
the  distances  from  the  dipoles  to  the  shields  are  d^  *  d^  ■  12",  as  shown  in 
Fig.  5. 

The  variations  of  the  shielding  coefficient  with  frequency  for  magnetic 
dipoles,  electric  dipoles,  and  plane  waves  are  shown  in  Fig.  6.  Experimental 
results  for  magnetic  dipoles  are  also  given.  From  these  results,  we  know  that 
the  copper  shield,  which  has  a  very  high  electric  conductivity,  possesses  a 
higher  shielding  coefficient  for  electric  dipoles. 

It  is  worth  noting  that  in  the  above  calculation  we  have  assumed  that  the 
propagation  direction  of  each  of  the  components  radiated  by  the  dipole  is  always 
perpendicular  to  the  boundaries  in  the  double  copper  shield  (see  Fig.  7).  This 
is  because  the  electrical  conductivity  of  copper  is  very  high. 


S  (dB) 


Figure  6*  Theoretical  and  experimental  results  for  electric  dipoles,  plane 
wave  and  magnetic  dipoles. 


6#  CONCLUSION 


The  formulas  given  in  References  [3]  and  [4]  are  only  for  specially 
oriented  dipoles,  i.e.,  the  dipoles  must  be  perpendicular  to  the  boundary.  Th 
formulas  are  not  available  for  arbitrarily  oriented  dipoles.  We  have  already 
obtained  formulas  for  z-oriented  and  x-oriented  dipoles.  The  former  is  perpen 
dicular  to  the  boundary,  whereas  the  latter  is  parallel.  If  a  dipole  is 
oriented  in  an  arbitrary  direction  to  the  boundary,  we  can  always  resolve  the 
dipole  into  two  equivalent  dipoles.  One  is  z-oriented,  and  the  other  is  x- 
oriented.  Therefore,  the  above  approach  can  be  used  for  arbitrarily  oriented 
dipoles.  Furthermore,  the  boundary  transformation  matrices  obtained  in  this 
paper  are  also  useful  for  the  study  of  wave  propagation  in  stratified  media. 
The  method  presented  in  this  paper  is  simple  and  convenient. 


APPENDIX  1 


FIELD  EQUATIONS 


Ue  have 


rA  +  kA  -  -Idi  5(x-x')  B(y-y')  <S(z-z') 


Let 


Az(x,y,z)  -  ///A  (a,S,y)  eJ  )+B(y-y'  )+y|z-z'  1 1  dad0dY 


Fourier  transforming  Eq.  (A.l)  we  get 


<k2  -  (a2  +  02  +  y2))  A  (a#6. 


Y) 


Id£ 


(2ir)- 


Therefore, 


Id* 


A2(o,B,y)  - - rr  2 - 2 - = - — 

(2*)3(k2  -  (a2  +  02  +  y2)) 


Thus  the  solution  of  Eq,  (A.l)  is 


*, (x.y.z)  ..Siifff  e3 1 «(*-*•  HB(y-y')]  doJsd , 


8ir3  - 


-  //  eJl°(x-x,)+e(y“y’)]  dod0  j 


*  ^jY(z-z') 


8*3  - 


(A.l) 


(A. 2) 


(k2  -  a2  -  B2)  -  Y2 


dY 


(A. 3) 


^et  k2  -  a2  -  B2  *  K2 ,  and  K  *  *^c2  ~  ot2  -  B2 ,  Then  the  second  integral  of 
Eq.  (A. 3)  becomes 


•  ejy(z-z') 

/ 


k2  -  y2 


(A. 4) 


Applying  the  complex  function  theory  to  evaluate  this  integration,  it  is  evident 

ejY(z-2’) 


that  the  Integrand  — - - -  has  two  poles  y  ■  on  the  real  axis  in 


the  complex  y-plane.  From  the  radiation  condition,  we  know  that  only  one  pole 


is  reasonable  and  that 


y  •  R(cos  $  +  jsin  $) 


Therefore, 


|eJ 2 * )  |  ■  e-R(z-z’)  sin  $ 


(A. 5) 


Since  sin  $  >  0  in  the  upper-half  plane,  the  integrand  will  become  negligibly 


small  as  R  increases  without  bound,  but  only  if  z  >  z'.  Thus,  we  can  take  the 


contour  as  shown  in  Fig.  A.l.  For  z  >  z',  the  contour  Includes  only  one  pole 


(r  ■  +K).  We  then  obtain  in  the  limit 


•  jy(z-z')  jr(z-z’) 

/  ■S“5 - 5“  dy  -  2irj  Res  -^—5 - 5—  y  ■  ±K 

—  K  -  y  K  -  y 


(A. 6) 


(z  >  z') 


1  2  ~2 
c  -  a  -  0 


Obviously,  for  z  <  z'  the  contour  must  be  closed  in  the  lower-half  plane. 


as  shown  in  Fig.  A. 2.  The  contour  includes  another  pole  (y  *  -K) ,  so  that 


•  Jy(z-z') 


/  m _ 

K2  -  y2 


-2  irj  Res 


-jr(z-z' ) 


K2  -  y2 


t2-a2-02(z-z' ) 


(A. 7) 


-  a2  -  02 


(z  <  z') 


and  for  z  >  z'  or  z  <  z' 


j  y( z-z  * ) 


J  'S“2 - 2“  dY 

—  *  -  yZ 


jA  -0  -3  | z-z* 


'l  2  2 

c  -  n  -  0 


z  >  z' 


or  z  <  z 


(A. 8) 


Finally,  we  have 


x  Id*  ?t  J 


j  [  a(  x-x f )+ 8(  y-y f )+ a2-  g2 


-a  -8  z-z' 


A2  -  a2  -  02 


dod8  (A. 9) 


v  ‘ v.v.v.v ■  v; 
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APPENDIX  2 


COORDINATE  TRANSFORMATION 


Considering  the  coordinate  transformation  from  the  E  ,  E  ,  E  coordinate 

x  y  z 

system  to  the  system  defined  by  the  orthogonal  triad  consisting  of  the  wave  vec¬ 
tor  k  and  vectors  E,  and  E,  shown  in  Fig.  A. 3  ,  we  obtain  a  coordinate  transfor- 

«  1 

mat ion  matrix 


/—  -n 

f-  -v 
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¥ 

X 

T. 

-  [R] 

¥ 
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J- 

where 


[R] 


sin  $ 


-cos  $ 


cos  $  cos  9  ,  sin  p  cos  9  ,  -sin  9 


The  final  results  are 
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(A. 10) 


(A. 11) 


(A. 12) 


(A. 13) 


(A. 14) 
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